Complex crystal structures formed by the self assembly of di-tethered nanospheres 



o 
o 

c 

a 

oo 



o 



I 

C 

o 
o 



> 
in 
p 

o 

On 
O 



X 



Christopher R. Iacovella 1 and Sharon C. Glotzer 1 ' 2 
1 Department of Chemical Engineering and 2 Department of Materials Science and Engineering 
University of Michigan, Ann Arbor, Michigan 48109-2136 
(Dated: January 8, 2009) 

We report the results from a computational study of the self-assembly of amphiphilic di-tethered 
nanospheres using molecular simulation. As a function of the interaction strength and direction- 
ality of the tether-tether interactions, we predict the formation of four highly ordered phases not 
previously reported for nanoparticle systems. We find a double diamond structure comprised of 
a zincblende (binary diamond) arrangement of spherical micelles with a complementary diamond 
network of nanoparticles (ZnS/D); a phase of alternating spherical micelles in a NaCl structure 
with a complementary simple cubic network of nanoparticles to form an overall crystal structure 
identical to that of AlCu2Mn (NaCl/SC); an alternating tetragonal ordered cylinder phase with a 
tetragonal mesh of nanoparticles described by the [8,8,4] Archimedean tiling (TC/T); and an alter- 
nating diamond phase in which both diamond networks are formed by the tethers (AD) within a 
nanoparticle matrix. We compare these structures with those observed in linear and star triblock 
copolymer systems. 



I. INTRODUCTION 

Recent attention in the literature has focused on meth- 
ods to self- assemble nanometer and micron sized par- 
ticles into highly ordered structures. For example, bi- 
nary nanoparticle superlattices [l[ and ionic colloidal 
crystals Q have been reported with structures reminis- 
cent of atomic crystals, e.g. NaCl and CsCl. Com- 
plex phases, such as the double diamond lattice ||, 
have been synthesized for covalently linked systems of 
tetrapods. Other methods to assemble nanoparticles 
have relied on anisotropic interactions, e.g. dip ole mo- 
ments @, H, H, 0] an d surface patterns [j| H, 03 EL E2] > 
to create wires0, H, EH , free-floating sheets U and 2-d 
crystals Q- In previous work, we examined two methods 
to assemble particles into ordered arrays and other com- 
plex structures using anisotropic interactions: "patchy 
particles" [12] and "tethered nanoparticles" [l4j. Patchy 
particles are nano or micron sized particles with direc- 
tional interactions conferred via patches on their surface. 
The location of patches dictates the local ordering and 
structure, ranging from square packed sheets [12| to the 
diamond lattice [la ]. Tethered nanoparticles are hybrid 
nanoparticle-polymer building blocks where nanoparti- 
cles are bonded as "head groups" to immiscible polymer 
tethers to create a new type of amphiphile. The immisci- 
bility between the nanoparticle and polymer tether facili- 
tates microphase separation into bulk periodic structures 
similar to those observed in block copolymers, including 
phases such as lamellar sheets 0, [l7], EaL ordered cylin- 
ders [HIISGJ], and the double gyroid~]ia,[l9j], but with 
additional ordering arising from the nanoparticle shape. 

Here we examine the self-assembly of di-tethered 
nanospheres (DTNS), nanospheres that have two short 
polymer tethers attached to their surfaces, separated by 
an angle 9 that we vary from 30° to 180°. We explore 
the phase behavior as a function of both immiscibility 
(via interaction strength) and directionality of the inter- 
actions (via the planar angle 9 between tethers). As a 



result of the immiscibility and directionality of the inter- 
actions, DTNS possess characteristics of both tethered 
nanoparticles and patchy particles. We explore a regime 
where tether length is short enough such that tether lo- 
cation is correlated to the attachment point, but not so 
short as to approach the limit of a patchy particle. A 
schematic is shown in figures QJa-c). 
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FIG. 1: Schematic of (a) patchy particle with localized 
patches on the surface, (b) DTNS with short tethers, and 
(c) DTNS with long tethers. The arrows are drawn from the 
center of the nanosphere through the attachment site to high- 
light the correlation between tether location and attachment. 
We conduct simulations in the range of schematic (b) where 
tether locations are correlated to their attachment points, (d) 
Schematic of the DTNS building block. The blue colored 
tether is labeled as A, the nanoparticle as B, and the red 
tether as C. 



II. MODEL AND METHODS 

We consider a general class of tethered nanoparticles 
rather than any one specific system and use empirical 
pair potentials that have been successful in the study of 
block copolymers, surfactants, and colloids, and in previ- 
ous studies of tethered and patchy nanoparticles. In the 
DTNS system, self-assembly is driven by immiscibility 
between the tethers and nanoparticles and controlled by 
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(1) the angle 8 between the tethers, (2) the size and shape 
of the nanoparticle, and (3) the interaction strength be- 
tween like species. We model the nanospheres as beads of 
diameter D=2.5cr connected to tethers via finitely exten- 
sible non-linear clastic (FENE) springs, where the max- 
imum allowable separation is set to Ro = 1.5a and the 
spring constant is set to k = 30. Tethers are modeled 
as bead-spring chains containing four beads of diame- 
ter a connected via FENE springs. The planar angle, 
8, between tethers at the surface of the nanoparticle is 
controlled by the use of a harmonic spring, with k = 30 
and Ro = (D + a)sin(6/2) . A schematic of the build- 
ing block is shown in figure [ljd) . To realize long time 
scales and large systems required to self-assemble com- 
plex mesophases, we use the method of Brownian dynam- 
ics (BD)|20]. Our simulations are performed under melt- 
like conditions where like species are attractive and unlike 
species are not attractive. To model the attractive inter- 
action between like species, we use the Lennard- Jones 
potential (LJ), which induces demixing below a certain 
critical temperature. The LJ potential is given by, 
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where e is the attractive well depth, U s hi ft is the energy 
at the cutoff, a is the parameter used to shift the interac- 
tion to the surface of the nanoparticles, and r cut = 2.5a. 
For interactions between tethers of like species we set a 
= and for nanoparticles we set a = 1.5a to properly 
account for excluded volume. Species of different type 
interact via a purely repulsive Weeks-Chandler- Andersen 
(WCA) soft-sphere potential to account for short-range, 
excluded volume interactions. The WCA potential can 
be described by the LJ equation (eqn. Q}, with U s hift = e 
and r cut = 2 1 / 6 er. For nanoparticle-tether interactions 
we set a = 0.75a and dislike tethers a — 0. Volume 
fraction, </>, is defined as the ratio of excluded volume 
of the beads to the system volume and the degree of 
immiscibility and solvent quality are determined by the 
reciprocal temperature, 1/T* = tjk^T. Further details 
of the model and method applied to tethered nanoparti- 
cles can be found in reference [lij ]. We emphasize that 
for this study, the potentials are chosen such that they 
capture the overall nanoparticle-tether and tether-tether 
immiscibility, the geometry of the nanoparticle, and the 
angle between the tethers. Changes to the phase behav- 
ior are expected if the individual interaction strengths are 
changed asymmetrically. Throughout the paper, figures 
are color coded as shown in figure [Ud) ; A tethers and 
the aggregates they form are colored blue, B nanoparti- 
cles are white/gray, and C tethers are red. 

We performed simulations at a volume fraction <j) = 
0.45 utilizing the following simulation procedure. For a 
given 8, we start from a high-temperature, disordered 
equilibrated state and incrementally cool the system, al- 
lowing the potential energy to equilibrate for several mil- 
lion time steps at each T* before cooling again. We sim- 



ulated over 150 individual statepoints to calculate the 
"phase diagram" presented below. For each 8 we per- 
formed at least two different cooling sequences to ascer- 
tain the path independence of the observed phases. We 
also performed simulations at different system sizes to 
eliminate artifacts due to box size. 

To identify the mesophases we utilize a combination of 
visual inspection, calculation of the structure factor [2lj |. 
and construction of the bond order diagram (BOD). The 
BOD shows the directions of all vectors drawn from a 
particle or micelle to neighboring particles/micelles pro- 
jected on the surface of a sphere, creating an "average" 
picture of the orientational order in the system. Sys- 
tems that have highly correlated neighbor directions (e.g. 
bulk crystalline materials) will show distinct groupings of 
points on the surface of the sphere; in contrast, a disor- 
dered system will appear as points randomly distributed 
on the surface. To analyze structures comprised of self- 
assembled micelles, we locate the center of mass of each 
micelle by modifying an image processing technique de- 
veloped by Crocker and Grier that is typically used to 
identify the center of colloidal particles from microscopy 
data [22fl . Substituting microscopy data by simulation 
data, we construct a density profile of our system and 
perform the standard Gaussian filtering and centroid cal- 
culation detailed in reference 12211 . 
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FIG. 2: Phase diagram for 1/T* vs. 9 showing four dis- 
tinct phases. From left to right: (1) a binary mixture of 
spherical micelles ordered into a zincblende diamond lattice 
of tethers with a complementary diamond network formed by 
the nanospheres (ZnS/D); (2) a binary mixture of spherical 
micelles of tethers ordered into an NaCl lattice with a com- 
plementary simple cubic network formed by nanospheres to 
form an overall crystal structure identical to that of AlCu2Mn 
(NaCl/SC); (3) a checkerboard pattern of alternating tetrag- 
onally ordered cylinders of tethers with a complementary 
tetragonal mesh of nanoparticles in an [8,8,4] Archimedean 
tiling (TC/T); and (4) an alternating diamond network (AD) 
of tethers within a nanoparticle matrix. Approximate phase 
boundaries are drawn to help guide the eye. Symbols refer to 
simulated state points; open blue squares refer to disordered 
states. 
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III. RESULTS 

We predict four phases not previously reported for 
tethered nanoparticle systems. These phases are indi- 
cated in the phase diagram shown in figure [2] and dis- 
cussed individually below. 




FIG. 3: (a) Centers of mass of the micelles formed by tethers 
that order into the ZnS lattice showing 27 unit cells with a 
unit cell size of approximately 8<j. A perfect diamond lattice 
is shown in the inset, (b-c) BOD for nearest neighbors in the 
ZnS structure for the two different tetrahedral arrangements, 
(d) The combination of figures (b) and (c). (e) Diamond net- 
work formed by the nanoparticles in the Zns/D phase, (f) 
S(q) for the nanoparticle strut in the binary diamond struc- 
ture (gray) for and S(q) for the perfect diamond structure 
(black). All DTNS data at 6 = 30° and 1/T* = 0.8. 



A. Zincblende ordered spherical micelles/diamond 
network (ZnS/D) 

For the range 30° < 9 < 45°, we find a binary mixture 
of spherical micelles ordered into a zincblende diamond 
lattice of tethers with a complementary diamond network 
formed by the nanospheres (ZnS/D). In figure 05 a) , we 
show a simulation snapshot of the centers of mass of the 
micelles formed by the tethers in the zincblende struc- 
ture. This structure is the two-component analog of the 
diamond lattice [23|. Figures [3^b)-(d) show the BODs of 
the micelles centers of mass for the zincblende structure. 
We split the BOD into two separate diagrams, since the 
diamond phase possesses two bond configurations that 
are 60° rotations of each other. We plotted clusters where 
a "blue" particle (i.e. micelle formed by the A portion) is 
at the center surrounded by "red" particles (i.e. micelles 
formed by the C portion), and a second diagram where 
clusters have a "red" particle at the center surrounded by 
"blue" particles; these cluster definitions properly group 
the data by orientation of the tetrahedrons. Both BODs 



in figure [3jb) and (c) show clear tetrahedral arrange- 
ments; the BODs for an ideal zincblende structure are 
plotted as lines in both plots, and agree with our simu- 
lation results. The blue centered and red centered tetra- 
hedrons have complementary orientations (i.e. rotations 
of 60°), as shown in figure [3jd) . 

Within this phase, the nanoparticles organize into a 
diamond network that is woven into the micellar lattice. 
An isosurface of the nanoparticle diamond network struc- 
ture is shown in figure[3]Je). In figure[3Jf) we plotted S(q) 
for the nanoparticles and a perfect diamond structure; 
the perfect structure was scaled such that the unit cell 
size matched our simulation. We find that the peaks cor- 
respond well to the perfect diamond structure; the small 
deviations can be attributed to thermal noise and the fact 
that our diamond network is composed of locally disor- 
dered nanoparticles (i.e. they have microphase separated 
into an ordered structure, but are not locally ordered). 
The overall ZnS/D phase is composed of two interwo- 
ven diamond structures where one network is formed by 
nanoparticles and the other consists of a binary lattice of 
spherical micelles. 



B. NaCl ordered spherical micelles/simple cubic 
network (NaCl/SC) 

Within the range 60° < 6 < 75°, we find a binary 
mixture of spherical micelles of tethers ordered into an 
NaCl lattice with a complementary simple cubic network 
formed by nanoparticles (NaCl/SC). In figure B|a) we 
plot a simulation snapshot of the centers of mass of the 
micelles formed by the tethers that order into a NaCl lat- 
tice. We can see that the structure clearly demonstrates 
alternating chemical specificity. Figure 0Jb) shows the 
BOD for the centers of mass of the micelles, where we ig- 
nore chemical specificity of the micelles and simply calcu- 
late the BOD for nearest neighbor micelles. The resulting 
BOD shows a simple cubic arrangement of the micelles, 
which corresponds to the BOD of a perfect NaCl system, 
shown as lines in[4jb). Figure IUc) shows the BOD for a 
perfect CsCl structure (often observed in triblock copoly- 
mers, as discussed later); it is clear that we have NaCl 
rather than CsCl. 

In this phase the nanoparticles fill the space between 
the micelles, aggregating at the interstitials and arrang- 
ing into a simple cubic network; an isosurface of the 
nanoparticles is shown in figure E^d) . The nodes of the 
simple cubic structure each have six connections points; 
a single node is shown in figure EJe) ■ The highest den- 
sity of nanoparticles (i.e. the node) resides in the center 
of an eight-particle NaCl unit cell, as shown in 0Je) - 
this high density location corresponds to the placement 
of the central particle in a BCC lattice. We utilized the 
same procedure to approximate the location of the nodes 
as we used to calculate center of mass of the spherical 
micelles. Figure 0|f) shows an eight-particle unit cell of 
NaCl extracted from our system with the locations of the 
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FIG. 4: (a) Centers of mass of the NaCl ordered micelles 
formed by the tethers; the unit cell size is approximately 6.5a. 
A perfect NaCl unit cell is inset, (b) BOD of the micelles cen- 
ter of mass; the BOD of a perfect NaCl structure is shown as 
lines, (c) BOD of a perfect CsCl structure for comparison, 
(d) Isosurface of the nanoparticles showing a simple cubic 
network arrangement, (e) 8 particle unit cell of NaCl formed 
by the micelles in the system, with the node of the nanopar- 
ticles network at the interstitial shown as an isosurface. (f) 8 
particle unit cell of NaCl formed by the micelles in the sys- 
tem with the nodes of the nanoparticle network drawn as gray 
spheres showing the AlCu2Mn structure. All data for 9 — 60° 
and 1/T*=0.8. 



cylinders. The nanoparticles organize into a 3-d structure 
whose cross section taken along the axis of the cylinders 
is a tetragonal mesh that separates the cylinders. The 
overall arrangement of species in the system can be de- 
scribed as the [8,8,4] Archimedean tiling, constructed of 
octagons and squares; this tiling is shown in figure E^b) , 
overlaid on the simulation data. 




FIG. 5: (a) End view of the TC phase at 9 = 90° and 1/T* 
= 0.8, where tethers form cylinders and nanoparticles form 
a square lattice that repeats along the long dimension, (b) 
[8,8,4] Archimedean tiling overlayed on tetragonal cylinders. 
All species are rendered at their true size, (c) Isosurfaces 
of the two tether domains, which self-assembled into the AD 
structure at 9 = 180° and 1/T* = 0.8 showing 27 unit cells 
with a unit cell size of 10<r; nanospheres are rendered at half 
their true size, (d) Structure factor of the AD phase. 



nanoparticle nodes rendered as gray spheres. The overall 
phase corresponds to the AlCu2Mn structure, also known 
as the Heusler (L2I) phase [23j. The AlCuaMn structure 
is a three-component analog of BCC (CsCl is the two- 
component analog to BCC). 



C. [8,8,4] Alternating tetragonal 
cylinders/tetragonal mesh (TC/T) 

For the range of 90° < 9 < 180° we find a cross- 
sectional checkerboard pattern of alternating tetrago- 
nally ordered cylinders of tethers with a complementary 
tetragonal mesh of nanoparticles to form an overall [8,8,4] 
Archimedean tiling (TC/T) (Fig. Ufa)). The cylinders 
alternate by the tether type that forms them, creating a 
checkerboard pattern when viewed along the axes of the 



D. Alternating diamond (AD) 

For the range of 150° < 9 < 180° we find an alternat- 
ing diamond network (AD) of tethers within a nanopar- 
ticle matrix. The AD phase, shown in figure [HJc) for 6 
= 180° and 1/T* = 0.8, consists of two chemically dis- 
tinct, interpenetrating diamond networks, one formed by 
the A tether and one formed by the C tether. Each di- 
amond network is composed of cylindrical tubes, where 
four tubes connect at a node in a tetrahedral arrange- 
ment. This structure was identified visually and by cal- 
culating the structure factor, plotted in figure^d), show- 
ing characteristic peaks with ratio y2 : \/3 : v4 : V6 as 
expected [24| : for ease of viewing the x-axis was scaled 
by 1.34, such that the numerical values on the x-axis cor- 
respond to the values in the characteristic ratio (i.e. the 
first peak occurs at a numerical value of V2) . 
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IV. DISCUSSION 

We can obtain insight into the observed DTNS phases 
by comparing with what is known about triblock copoly- 
mers. The DTNS system can be loosely thought of as 
a nanoparticle equivalent of an ABC triblock copolymer 
where, in this case, the center block of a triblock has been 
replaced by a nanosphere. As we decrease 9 between the 
two tethers in the DTNS system, we see a change from 
alternating diamond to alternating tetragonal cylinders 
to spherical micelle phases, as shown in figure [2] In lin- 
ear triblock copolymers, as the length of the middle block 
(B block) is decreased, we typically see a change from a 
tricontinuous structure (e.g. the alternating gyroid (2f| 
or alternating diamond phase [26|, [53, HH ) to alternating 
tetragonal cylinders [25|, [28| to CsCl structured micelles 
[25l . In both cases, the net effect is that the A and 
C blocks are brought closer together, constraining the 
possible tether configurations and resulting in a phase 
transition. In a rough sense, the overall progression of 
phases is similar whether we change the angle between 
tethers or, in linear triblocks, decrease the length of the 
middle block. 

The AD and TC/T phases are both well known in 
the linear triblock copolymer literature and their forma- 
tion in the DTNS system is not entirely surprising; in 
these two cases, the DTNS behave very similar to lin- 
ear triblocks and the geometry of the nanoparticle ap- 
pears to have little impact on the resulting structure. 
However, as we noted earlier, linear triblock copolymers 
have been shown to form CsCl ordered spherical micelles 
rather than the NaCl ordered micelles we find for DTNS. 
Furthermore, it has been calculated that for linear tri- 
block copolymers CsCl has a lower free energy than NaCl 
[29l | . thus, we would expect that DTNS might also form 
the CsCl structure. In previous work |16j, we found 
that under selective solvent conditions, mono-tethcrcd 
nanospheres could be roughly mapped to surfactants and 
block copolymers by considering the ratio of the excluded 
volume of head to tail, [161 ], however, even when a 
nanosphere and polymer chain have the same excluded 
volume, there can be important, subtle differences. A 
nanosphere has fewer configurational degrees of freedom 
than an equivalent polymer with the same excluded vol- 
ume; the radius and effective volume (i.e. shape or mass 
distribution) of the nanosphere are constant, whereas the 
radius of gyration and effective volume of a flexible poly- 
mer can vary based on solvent conditions, temperature 
and volume fraction. As a result, the flexible middle 
block in a triblock copolymer may reduce the correla- 
tion between the A and C blocks, e.g. by bending or 
stretching, allowing the A and C blocks a larger config- 
uration space. In the DTNS system, we remove many of 
the degrees of freedom of the middle block by replacing 
it with a nanosphere that has a fixed volume/geometric 
contribution and by including a bond angle constraint, 
both of which limit the configurational entropy of the 
A and C tethers. As such, it is reasonable to conclude 



that for the DTNS system, these changes may result in a 
different favorable configuration of the individual blocks, 
thus changing the overall global structure, in this case 
stabilizing the NaCl structure over CsCl. 

The effect of architecture of the DTNS (via the pla- 
nar angle) is clear when we examine the transition from 
NaCl/SC to the ZnS/D. As the spacing between the first 
beads of the two tethers approaches ^lcr, the tethers 
are essentially connected to the same location on the 
nanosphere, increasing the prominence of the nanosphere 
and making it less linear with respect to the tethers; in 
this limit the DTNS more closely resembles a star tri- 
block copolymer with one collapsed block, as sketched 
in figures [^a-b). This change in architecture brings 
the tethers closer together which in turn brings the mi- 
celles closer together; the average distance between the 
centers-of-mass of the micelles in the NaCl configuration 
is 6.6I0- ± 0.57 for 9 = 60° and 1/T* = 0.8, while in 
the zincblende configuration the distance is 5.86<r ± 0.91 
for 9 = 30° and 1/T* = 0.8 (note: both data sets are 
Gaussian and the size of the micelles are equivalent be- 
tween the two systems). By bringing the micelles closer 
together we force a change in the aggregation behavior of 
the nanospheres, which can be observed not only be the 
change from a simple cubic network to diamond network, 
but also by the difference in nanoparticle coordination 
where cn = 6.86 ± 1.24 for NaCl/SC verses cn = 7.46 ± 
1.37 for Zns/D. Star triblock copolymers are not known 
to form the ZnS/D structure, but are known to exhibit 
a columnar phase described by the [6,6,6] Archimedean 
tiling [30(. In the [6,6,6] phase, all three blocks form 
cylinders that are arranged in an alternating hexagonal 
pattern (30j . as shown in figure [S^c). If we only consider 
two of the three types of hexagons (e.g. red and blue), 
essentially excluding the portion that would correspond 
to the nanoparticles in our system, we find a binary ar- 
rangement of hexagons in six-member rings that closely 
resemble the structure of graphite. This [6,6,6] structure 
is nearly identical to the 2-d projection of the ZnS/D 
structure, as shown in figure Od) where three graphite 
hexagons are highlighted. In other words, the ZnS/D 
structure formed by the spherical micelles is essentially 
the 3-d micellar equivalent of the columnar [6,6,6] phase. 
Just as in the NaCl/SC phase, the geometry and volume 
contribution of the nanoparticle seems to be important 
in stabilizing the ZnS /D structure over the typically ob- 
served triblock copolymer phase; the bulky geometry of 
the nanoparticle helps to induce curvature, stabilizing 
spherical rather than cylindrical micelles. If we decrease 
the diameter of the nanoparticle from D=2.5tr to 2. Oct, we 
reduce the ability of the nanoparticle to induce curvature. 
We simulated a system with D=2.0cr at 9 = 30° finding 
the [6,6,6] columnar phase, as shown in figure[B{e) rather 
than ZnS/D. If the length of the B block for a star tri- 
block copolymer in the [6,6,6] phase were increased, simi- 
lar to changing the nanosphere diameter from D=2.0cr to 
2.5a, the system has been shown to form the [8,6,4;8,6,6] 
columnar phase [3(|, not ZnS/D. This clearly highlights 
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the importance of the geometry of the nanoparticle; it is 
not only the volume of the nanoparticle that is important 
but also the distribution of volume (i.e. the shape). 




FIG. 6: (a) Schematic of a star triblock copolymer, (b) 
Schematic of the DTNS building block at 0=30°. (c) 
Schematic of the [6,6,6] Archimedean tiling. A graphite struc- 
ture is shown as the black hexagonal lattice, (d) Three 
equivalent graphite tiles are overlaid on the projection on the 
zincblende structure, (e) Reducing the size of the nanoparti- 
cle from D=2.5ct to 2. Oct results in the formation of the [6,6,6] 
columnar structure. Cylinders of tethers are shown as isosur- 
faces for clarity. 



V. CONCLUSIONS 



Our results show that the structural phase behavior of 
di-tethered nanospheres is a function of both the direc- 
tionality and strength of the tether-tether interactions. 
We have shown that DTNS can produce unique struc- 
tural arrangements of both tethers and nanoparticles and 
that these arrangements are controlled by the geometry 
of the nanoparticle and location of the attached teth- 
ers. We have demonstrated a novel route to form dia- 
mond and SC networks of nanoparticles, two structures 
highly sought for photonics applications [31| • Overall, we 
have shown that the use of soft-matter tethers with direc- 
tionality can be used to produce highly ordered periodic 
structures that would not necessarily be expected of ei- 
ther equivalent flexible polymer systems or pure nanopar- 
ticle systems in the absence of tethers. 
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